Measuring roughness of buried interfaces by sputter depth profiling 
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In this communication, we report results of a high resolution sputter depth profiling analysis of 
a stack of 16 alternating MgO and ZnO nanolayers grown by atomic layer deposition (ALD) with 
thickness of ~5.5 nm per layer. We used an improved dual beam approach featuring a low energy 
normally incident direct current sputtering/milling ion beam (first beam). Intensities of 24 Mg + and 
64 Zn + secondary ions generated by a pulsed analysis ion beam (second beam) were measured as a 
function of sample depth by time-of-flight secondary ion mass spectrometry (TOF SIMS). 

Experimental results of this dual beam TOF SIMS depth profiling processed in the framework 
of the mixing-roughness-information (MRI) model formalism demonstrate that such an approach is 
capable of providing structural information for layers just a few nm thick. Namely, it was established 
that the interfacial roughness of the MgO/ZnO multilayer structure equals 1.5 nm. This finding 
by TOF SIMS was cross-validated by independent measurements with specular X-ray reflectivity 
(XRR) technique. In addition, the TOF SIMS-MRI analysis suggests that the obtained 1.5 nm 
roughness should be attributed to the native roughness (jagged type) of the interface rather than 
to interdiffusion at the interface during the ALD synthesis. 

PACS numbers: 68.35.bg, 68.35.Ct, 81.15.Gh, 82.80.Ms, 82.80.Rt 



I. INTRODUCTION 

Quality of interfaces in multilayered materials is 
of great interest and importance, since the way lay- 
ers adhere to each other greatly affects their perfor- 
mance: spin-based electronics^ giant magnetoresistance 
structures/devices? 2 - high-fc multilayer materials^ artifi- 
cial Bragg crystals,— heat transfer in nanostructures,— and 
others. 

A powerful and efficient set of methods commonly used 
to determine (i) single layer thickness, (ii) roughness and 
(iii) density in a multilayer system combines the X-ray 
reflectivity (XRR) method with ellipsometry. XRR tech- 
nique uses the values of (i), (ii), and (iii) as fitting pa- 
rameters jointly with some type of boundary conditions 
to solve a general inverse problem of wave scattering 
XRR results can be partially validated by optical ellip- 
sometry, which is also based on solving a similar inverse 
problem. It relies on refractive index (a function of den- 
sity) and a layer thickness in a single/multiple layered 
structure, and depends on boundary conditions as well. 
This makes the data fitting model sensitive and depen- 
dent on many free parameters, which can grow fast when 
the number of layers in the layered structure increases. 
Moreover, such indirect techniques based on measuring 
optical responses of a material may become inefficient if 
one needs to work with three-layer stacks (or even more 
general case of n layers) where repetition of certain lay- 
ers in an n-layer stack varies, as, for instance, occurs 
in the case of metamaterialsi^ A problem will also ex- 
ist for multilayer semiconductor heterostructures, where 
the composition variation between nearest layers is about 
the same with a few percent difference in solid solution 
content,- Hence, the change in density is so small that 
the sensitivity of the XRR method may be insufficient for 
characterization of such materials. Thus, the search for a 



more direct and established approach to characterization 
of nm-thick layers in a variety of multilayered structures, 
which is sensitive and independent on optical properties, 
is highly desirable. It can complement XRR and ellip- 
sometry techniques, and in some cases might do the job, 
which is beyond their capabilities. 

Mass spectrometry retains its importance in basic 
physics fii chemistryfi 2 . and in materials science 1 ^ due to 
a number of unique advantages, such as sensitivity to iso- 
topic ratio, trace amounts of elemental/molecular species 
on the surface and in the bulk, etc. In this context, an 
extension of its capabilities is of great interest. In the 
present paper, we demonstrate that properly designed 
mass spectrometry of secondary ion species in sputter 
depth profiling experiments can successfully reveal struc- 
tural features of nanolayered structures and characterize 
their surface and interfaces. To this end, an improved 
dual beam secondary ion mass spectrometry (SIMS) with 
low energy normal incidence sputtering/milling was ap- 
plied to depth profile a layered MgO/ZnOx8 structure 
grown by atomic layer deposition. The obtained depth 
profiles were processed in terms of mixing-roughness- 
information (MRI) model. This model has analytical so- 
lutions with, in case of SIMS, two variable parameters 
having straightforward physical interpretation: the ion 
beam mixing length and the roughness. The SIMS-MRI 
results are compared to structural parameters obtained 
by specular XRR measurements to confirm the correct- 
ness of made conclusions. 



II. SAMPLES AND EXPERIMENTAL 

A layered structure |5.5-nm MgO/5.5-nm ZnO| x8 was 
grown by ALD on a Si substrate using established pre- 
cursor chemistries for MgO 1 ^ and Zn0 15 i 16 and charac- 



2 



terized by X-ray diffraction and ellipsometry. The struc- 
tural data obtained from these measurements indicated 
that MgO layers are amorphous and ZnO layers are poly- 
crystalline in the wurtzite phase. The roughness of the Si 
substrate was ~0.3 nm. The initial layer-to-layer mixing 
due to thermal diffusion during growth at T=473 K is ex- 
pected to be extremely low. Specular X-ray reflectivity 
measurements were made with a Philips X'Pert Pro MRD 
diffractometer using Cu Ka radiation (A=1.5418 A) and 
operated at 30 kV/40 mA. The incident X-ray beam was 
conditioned by a 60 mm graded parabolic W/Si mir- 
ror with a 0.8° acceptance angle and a 1/32° divergence 
slit. The reflected beam was collected with a PW3011/20 
sealed proportional point detector positioned behind a 
0.27° parallel plate collimator. 

Secondary ion mass spectrometry studies were per- 
formed in a custom-designed SAPJSA (Surface Analysis 
by Resonance Ionization of Sputtered Atoms) instrument 
in the Materials Science Division at Argonne National 
Laboratory^ SARISA combines two independently op- 
timized Ar + ion beams: one of low energy (a few hundred 
eV) and normal incidence for ultimate depth resolution 
ion milling/sputtering, and another for elemental time- 
of-flight (TOF) SIMS analysis with high lateral resolu- 
tion. This arrangement is based on the known power- 
ful dual beam approach to depth profiling^ and can be 
dubbed gentleDB . There are three main advantages of 
this setup. First, depth resolution is controlled by the 
milling beam if the parameter a is much smaller than 
unity, or a -C 1 (see below). Second and third, normal 
incidence permits varying the impact energy of primary 
ions by target bias and does not introduce additional 
roughening of a sample surface! 19 i 20 

Elemental depth profiles, concentration/intensity ver- 
sus depth, of the sandwich structure were obtained by 
a sequence of alternating cycles. Ion milling by a raster 
scanned primary direct current Ar + ion beam at 500 eV 
and normal incidence^ 7 , is followed by TOF SIMS anal- 
ysis of revealed subsurface (various depths, from surface 
down to substrate) by a raster scanned pulsed (200 ns 
long) Ar + ion beam at 5 keV energy and 60° incidence 
with respect to the target normal, as shown schemati- 
cally in Fig. la. This pair of cycles is repeated multiple 
times until the Si substrate is reached, which is moni- 
tored by Si + peak intensity in a mass spectrum. In the 
depth profile, the Si substrate grows in as a spike at- 
tributed to ^100-fold enhanced secondary ion yield due 
to presence of Si02- It then stabilizes at a lower con- 
stant level, as shown in Fig. lb. The depth resolution is 
controlled ultimately by the milling beam characteristics. 
This condition is fulfilled if the parameter called effective 
erosion efficiency E = Y ■ j ■ t/e specific for the milling 
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Figure 1: (a) Simplified experimental setup: 8 MgO 
(cyan)/ZnO (red) pairs, superimposed and centered blue and 
green squares are craters from a raster scanned direct cur- 
rent normal incidence milling Ar + beam (impinging solid blue 
line), and a raster scanned area of a pulsed probing Ar + beam 
(impinging green dotted line), respectively, (b) A full gen- 
iZeDB SIMS depth profiles of the MgO/ZnOx8 ALD struc- 
ture on a Si substrate obtained at the 500 eV milling Ar + 
beam combined with the 5 keV probing Ar + beam. 



beam is much higher than that of the analytical beam, or 
a = E a /E s <C 1. Here, Y is a sputtering yield, which de- 
pends on a beam energy s, the primary projectile species 
and an angle they attack the surface, j is the current 
density, t is the total sputtering time during which an 
ion beam is on. In our experimental setup a was ~10 -5 . 

The milling beam was digitally raster scanned over a 
square area of ~1 mm 2 (blue square in Fig. la), while 
the analytical beam was raster scanned over a square of 
^500x500 /im 2 (green square in Fig. la). Both raster 
areas were precisely overlapped by using an in situ 
Schwarzschild microscope, 17 and ex situ by white light 
profilometryi^i 

Ion beam currents were measured and focused in situ 
by a custom made graphite Faraday cup (FC) consisting 
of an internal pin and inlet holes of 250 /jm dia. on the 
mask. Such a FC design permits to control focusing con- 
ditions of an ion beam. Ion beam profiles were checked 
and known to have a symmetric Gaussian distribution by 
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burning dents in soft materials and profiling the dents by 
ex situ white light profilometry^i It allows one to pre- 
cisely calculate a current density of an ion beam using a 
known dc current value measured on the FC. The milling 
beam parameters at 500 eV were 1 /iA dc and a FWHM 
of Gaussian-like beam profile of ~150 /im. The analyt- 
ical beam characteristics at 5 keV were 300 nA dc and 
a FWHM of Gaussian-like beam profile of ~40 fim, and 
then reduced down to 30nA by choosing an appropriate 
aperture. 

The measured energy spread Ae of our low energy 
column^ is 23 eV yielding Ae/e ~ 5 x 10~ 2 at 500 eV 
milling energy. Thus, the milling beam can be considered 
as monoenergetic with high precision. 



III. RESULTS, DISCUSSION AND 
CONCLUSIONS 

Fig. lb demonstrates SIMS depth profiles of Mg + and 
Zn + obtained by gentleDJi at 500 eV ion milling, re- 
flecting the full set of peaks in the periodical structure, 8 
peaks for MgO and ZnO each. Fig. 2a is a high-resolution 
profile obtained under the same conditions as the profiles 
in Fig. lb, but at small ion milling increments, so that the 
depth difference between two consecutive points is ^0.2 
nm, corresponding to monolayer thickness for wurtzite 
ZnO. 

We used the high-resolution Zn + depth profile for fur- 
ther evaluation of structure features of the multilayer 
ALD stack. This was accomplished in the terms of an 
analytical approach called mixing-roughness-information 
(MRI) modeh 22 i 23 In essence, the model has 3 physically 
meaningful parameters. Mixing is characterized by the 
ion beam mixing length w, a scale at which two per- 
fect layers of different compositions with perfect abrupt 
interface experience mutually uniform interpenetration 
under nonreactive gas bombardment. Roughness is the 
root-mean-square (rms) roughness a. By the term "infor- 
mation", the depth of collecting information A is meant, 
which is the same as the escape depth commonly used 
in surface analysis techniques. In general, A is a vari- 
able parameter, but in the case of SIMS it can be fixed 
at a value of 1-2 monolayers^ The effects of the two 
main processes, mixing and roughness, on the resulting 
peak depth profile are shown in Fig. 2b. It is seen that, 
in the case of an ideal rectangular layer, the ion beam 
induced mixing leads to asymmetric depth profile, which 
may get additional broadening and symmetrization due 
to the rms roughness, depending on the interplay between 
w and a values. The effect of the information depth A is 
not incorporated, since it is expected to be negligible. 

The green solid curve in Fig. 2a superimposed on the 
Zn + depth profile is the best fit of the experimental curve 
by the MRI model. The fit parameters were as follows: 
w is 0.4 nm and a is 1.5 nm, while A was fixed at 0.2 nm, 
corresponding to the wurtzite ZnO monolayer thickness. 
The flat thickness of the layer d, which is a thickness of an 
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Figure 2: (a) High-resolution gentleDH SIMS depth profiles 
of MgO/ZnO top layers at 500 eV milling Ar + beam and 5 
keV analytical Ar + beam. Mg signal is shown in cyan (open 
diamonds), and Zn is shown in red (open circles). The green 
solid curve is the best MRI model fit. (b) Principles of the 
MRI model describing how different distortion sources affect 
an ideal, rectangular (flat) layer profile of thickness d. The 
effect of the escape depth parameter, negligible in the case 
of SIMS, would slightly broaden the red curve (open squares) 
and is not shown. 



ideal flat layer (see Fig. 2b) defining boundary conditions 
for the model, appeared to be 3 nm yielding the total 
single ALD layer thickness to be of 6 nm. The sharpness 
of the last interface between the Si substrate and the ZnO 
layer can also be estimated. An up-slope of the Si signal 
shown in Fig. lb results in a roughness of ~0.4 nm. 

To cross-validate the results obtained by MRI mod- 
eling of the experimental data, we performed specular 
XRR measurements on a MgO/ZnO sample prepared un- 
der identical conditions. The use of XRR in this case 
seems to be a valid approach, since its sensitivity to fluc- 
tuations of the thickness (that can range in character- 
istic length scale) at the wavelength of ~1.5 A can be 
anticipated to be similar to that of SIMS based on a 
single ion penetration in the fine interfacial structure. 
Fig. 3 represents a typical specular XRR profile obtained 
for the MgO/ZnOx8 sandwich. Fitting of experimental 
data was performed with a commercial software pack- 
age (Panalytical X'Pert Reflectivity) which makes use of 
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Figure 3: A specular XRR profile of the MgO/ZnOx8 ALD 
structure on the Si substrate obtained at 1.5418 A irradiation. 
The black circles are the experimental reflectivity curve, while 
the simulated reflectivity curve is represented by the red solid 
line. 

the Parratt recursion formalism for reflectivity^ Due to 
the relatively large number of discreet layers present in 
the MgO/ZnOx8 heterostructure, several periodic con- 
straints were applied to reduce the set of fitting param- 
eters which would otherwise be required. A simple ide- 
ally periodic model assuming identical ZnO/MgO bilay- 
ers failed to reasonably approximate the data. A sec- 
ond model where the first ZnO layer and last MgO layer 
were treated discreetly as described by Jensen et al£ also 
failed to produce adequate fits to the experimental data. 
However, simulations where the density, p, thickness, and 
roughness, a, of each subsequent ZnO and MgO layer 
were constrained to increase linearly with increasing dis- 
tance from the substrate provided a good match to the 
experimental reflectivity curve. The experimental fit was 
further optimized by allowing each p, thickness, and a to 
vary from the linear model within reasonable physical 
constraints (e.g. bulk density values). The resulting re- 
flectivity fits confirmed the parameters values extracted 
from SIMS data: the average thicknesses of ZnO and 
MgO appeared to be 6.4 nm and 6.1 nm, with average 
roughness values of 1.1±0.2 nm and 1.6±0.3 nm, respec- 
tively. Parameters of the Si02 layer on the top of the Si 
substrate were also of the expected values: thickness of 
1.2 nm, <7=0.3 nm. 

Thus, XRR data independently proves that SIMS anal- 
ysis along with the chemical information is capable of re- 
vealing detailed internal interface structure, when distor- 
tions/artifacts related to a profiling procedure are min- 
imized to a scale, much smaller than that of a feature 
to be investigated. In the case of gentleDH approach, 
this means that the extracted mixing parameter w=0A 
nm is of the same order as the inherent escape depth of 
a secondary ion, a physical limitation of the technique. 



Therefore, the found 1.5 nm interfacial roughness of the 
MgO/ZnO layers is due to an ultra-short ion beam mix- 
ing length/lattice disturbance in gentleDH profiling pro- 
cedure. XRR also proves that the obtained 1.5 nm rough- 
ness is reasonable and falls well within the range provided 
by other experimental methods, both in this study and 
other studies of layered laminate oxides performed by 
AFM/STM^ XRR&Z and TEM4 

Let us extend the finding by gentleDB SIMS in con- 
junction with MRI results. First of all, in the context of 
the advantages of the low energy orthogonal sputtering 
described above, the rms roughness a is believed to be a 
native interface roughness of the grown ALD structure. 
In addition, it is known that XRR cannot distinguish be- 
tween the interfacial roughness and the interdiffusion. 7 
One may assume that the interdiffusion at the interface 
is a process similar in some way to the ion beam mix- 
ing, which normally causes a mutually uniform layers 
interpenetration, 23 thus blurring the initially sharp in- 
terfaces and forming a mixed-up interlayer. If the in- 
terdiffusion is pronounced, the MRI model suggests that 
the gentleDB SIMS depth profiles (keeping in mind the 
limit of w going down to zero) would look more asym- 
metric, closer to the green curve shown in Fig. 2b. For 
this particular MgO/ ZnO structure grown at only 473 K, 
the shapes of the Mg + and Zn + peaks, which are broad 
and nearly symmetric, in measured depth profiles sug- 
gest that the estimated 1.5 nm roughness corresponds to 
native roughness (i.e., "jagged" junction), rather than to 
interdiffusion, or that the interdiffusion does not exceed 
0.4 nm from the total 1.5 nm. 

To summarize, gentleDH SIMS, in conjunction with 
the MRI model, is a direct, but destructive and rela- 
tively slow method, which does not require prior assump- 
tions regarding the sample. On the contrary, the specular 
XRR is indirect, but nondestructive and relatively fast 
method, which is model dependent and needs prior infor- 
mation about the sample. Thus, when both techniques 
are applied to characterize buried interfaces in nanolay- 
ered materials, they complement and enhance each other, 
and yield comprehensive physical/structural and chem- 
ical information on the material: interfacial roughness 
distinguishing interdiffusion from native /jagged rough- 
ness, thickness, density, as well as elemental depth pro- 
files, including trace impurities, redistribution of isotopes 
in isotopically modulated structures^ and many other 
possible applications. 
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